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Using a relativistic Dirac-Brueckner analysis the OPEP contribution to the 
ground state energy of nuclear matter is studied. In the study the pion is derivative- 
coupled. We find that the role of the tensor force in the saturation mechanism 
is substantially reduced compared to its dominant role in a usual nonrelativis- 
tic treatment. We show that the damping of derivative-coupled OPEP is actu- 
ally due to the decrease of M* /M with increasing density. We point out that 
if derivative-coupled OPEP is the preferred form of nuclear eff'ective lagrangian 
nonrelativistic treatment of nuclear matter is in trouble. Lacking the notion of 
M* it cannot replicate the damping. We suggest an examination of the feasibility 
of using pseudoscalar coupled ttN interaction before reaching a final conclusion 
about nonrelativistic treatment of nuclear matter. 

The talk is based on a recent paper written by Tjon and myself [Q. 
The one pion exchange potential, OPEP, in momentum space is given by the 
expression: 

{\P + P'. \P-P'\^\\P + P. \P-P) = 



ff^AfAfm^ 2 (Ti ■ {p' ~ p)a2 ■ jp' -p) _ 



With ^( ^"^j^/"' )^ ~ llMeV the OPEP in coordinate space is given by the 
expression: 
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where the last term is the tensor force and 

<S'l2 = 3(Ti • f(T2 • f - (Ti • (72 

is the tensor operator. In a Hartree-Fock calculation in nuclear matter the S{r) 
potential is capable of contributing at saturation density, po — ^ni^, as much as ~ 
20MeV. But we must ignore it as it will be wiped out by the short-range correlation. 
The Yukawa potential and the tensor force contribute a mere ~ — 2MeV. Really 
important contributions of the tensor force come from second and higher orders. 
The large matrix element, {^Di \ S12 P Si) = VS, shown in the matrix below, 
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is responsible for this feature. 

The effect of the tensor force and its dominance in nonrelativistic nuclear physics 
are seen most dramatically from the following results for the deuteron 

{Deuteron \ T^cntiai | Deuteron) ^ 

I T4ensor P ^l) ~ -22MeV 

{Deuteron \ \ Deuteron) ^ +20MeV 

In a nonrelativistic Bethe-Brueckner calculation of nuclear matter one finds 
typically |l 

(N Matter | K I N Matter) „o„rciativistic ^ ~34 (p/po)"'*' MeV. (1) 

The exponent of the density p is markedly less than the nominally expected value 
of 1 because of Pauli blocking. 

Relativistic nuclear physics is heavily based on summing Bethe-Salpeter ladders 
using some form of Blankenbecler-Sugar-Logunov - Tavkhelidzc ^ prescription to 
obtain quasipotentials. In general quasipotentials are not simply described by tree 
graphs with dressed vertices. The OPEP is an exception. Other terms have the 
characteristics of four-point functions. In practice, one uses OBEP forms with form 
factors for approximate representation of quasipotentials. The parameters are fixed 
by fitting NN data. 

The relativistic results for the contribution of V-,^ to the deuteron is 

(Deuteron | V^, \ Deuteron) relativistic = -22MeV, 

suggesting an equally important role of the pion. In sharp contrast, this seems not 
to be the case in a relativistic treatment of nuclear matter. 



banerjee printed on February 9, 2008 



3 



Strong scalar {S) and vector (V_) fields of the order of a few hundred MeV are 
typical for relativistic theories 0, D based on a meson theoretical description 
of the nuclear force. These values are consistent with expectations based on the 
studies of scattering of ~ 1 GeV protons by nuclei. The large scalar fields have 
far reaching consequences in nuclear matter through the strongly medium modified 
nucleon mass M* = M + S. The saturation mechanism is believed to rest upon 
the decrease of magnitude of S with increasing density. Of course, in a mean field 
theory (MFT) like the QHD it is the only possible mechanism for saturation. 

The contributions of a particular meson exchange potential, Va , to the binding 
energy can be calculated using the Hellmann-Feynman theorem 

(N Matter | K | N Matter) ^ gl^^ ^ 2 (^M) • (2) 

"9aNN 

We find the following for the pion field contribution to E/A: 

(N Matter | \ NMatter)Reiativi.tic ~ -20 (p/po)°-'^ MeV. (3) 

^From this we see that the pion contribution is considerably suppressed compared 
to the value given by Eq. (||) for the nonrelativistic case. We will make clear that 
the suppression of OPEP is generic and not particular to the present calculation. 
Furthermore, OPEP has only a minor role in the saturation mechanism. This 
is exhibited in Fig. |l| where we plot our calculated results E/A (curve a) and 
E/A - {NM I I NM) - 17 MeV (curve b). The two curves have practically the 
same density dependence verifying that OPEP contributes little to the saturation 
mechanism. The subtraction of 17 MeV in curve (b) makes the scale more compact. 
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Fig. 1. Plots of tlie Dirac-Brueckner predictions of E/A (curve a) and 
E/A-{NM 1 I 7VM>-17MeV (curve b). 

The preceding numbers confirm the belief that in a relativistic treatment of 
nuclear matter the tensor force does not have the dominant role that it has in 



4 



BANERJEE PRINTED ON FEBRUARY 9, 2008 



the usual nonrelativistic treatment. The radicaUy different explanations of the 
saturation mechanism in nonrelativistic and relativistic studies of nuclear matter 
constitute a puzzling issue. A valid nonrelativistic treatment must reproduce the 
main physics of a valid relativistic treatment in leading order in w/c. Although the 
issue is a longstanding one, no resolution of it has been given to date. We have 
addressed this question. A Dirac-Brueckner (D-B) analysis || is at present the 
best tool we have for a relativistic study of nuclear matter. We examine here the 
role of derivative-coupled OPEP in D-B and show that it is substantially reduced 
due to relativity. Since the contribution of OPEP to the deuteron binding energy 
remains large in a relativistic treatment the damping in nuclear matter must be 
due to many-body effects. We find that it can be attributed to the decrease of 
M* /M with increasing density. 

The above results can be understood qualitatively by examining the second 
order contributions to the G-matrix. Keeping only the positive energy M* state 
contributions in the intermediate states we have, 
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x(p' I V{P) I fc)-i {k\V{P) |p).(4) 

The tensor force contributes mainly to the second term of Eq. (^) which makes 
the structure of the two-nucleon propagator important. Normally the boson masses 
in OBEP provide the scales for momenta in nuclear physics. But here we notice 
that the two-nucleon propagator provides a new scale, viz., M* . To exploit this 
new scale let us introduce dimensionless momenta, £ = p/A4*, n = k/M*, etc., 
and exhibit a few OBEP matrix elements in terms of these. 

a Exchange 



{■P\va\ k) 



[p— ky + ml 



M*2 (f_,-^)2 + (to^/M*)2 



(5) 
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Table 1. Table of AI* factors for various OBEP potentials. The number of stars in the 
2nd column indicates the importance of the OBEP in nuclear interaction. 



Parity 


Importance 


Boson 


Scaling Factor 


E 


^ ^ ^ 


(7 


(1/A/*)^ 


V 


^ ^ ^ 


Lu Dirac 


(1/Af*)2 


E 




p Dirac 


(1/Af*)^ 


N 




S 


(1/M*)2 





* * 


TT Der. Coupled 


(M*/M)2(l/M*)2 


D 


* 


p Pauli 


(M*/M)2(l/M*)2 


D 




rj Der. Coupled 


(M*/M)^(l/M*)^ 






LU Pauli 


(M*/M)2(l/M*)2 



TT Exchange (Derivative Coupling) 



dc , A ^ / 9^ ^ CTi • (A: - p)a2 ■ (k - p) 



2M' - (p*_/t)2+^2 



mt 



TT Exchange (Pseudoscalar Coupling) 
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r(Tr) ^1 ■^2-^^— — - — — — . (7) 



Notice the differences in the M* factors. The derivative-coupled pion-exchange 
potential has an extra damping factor of {M* /MY relative to the sigma-exchange 
potential. It is reasonable to expect that the M* /M factor suppresses derivative- 
coupled OPEP, Pauli coupled p, etc. 

No such damping factor is present for the pseudoscalar coupled pion-exchange 
potential. To our knowledge all published relativistic nuclear matter calculations 
have used derivative-coupling. Of course, the reason is well-known. The pair 
suppression problem is automatically taken care of with use of derivative-coupling. 
We will proceed as if derivative- coupling is correct. A discussion of pseudoscalar 
coupling vs. derivative-coupling follows at the end of the talk. 

The M* /M suppression is corroborated in more detail by the following calcu- 
lation. Let us modify the S obtained from the self-consistent D-B calculation by 
multiplying it with the factor a <1 thus generating a Al* = M + aS. By using the 
modified scalar self-energy in the nucleon propagators we recalculate first the G 
matrices and then E/A and finally (N Matter | | N Matter) using Eq. (1). Only 
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the a = 1 analysis is self consistent; others are not. But such a calculation is par- 
ticularly suitable to exhibit the role of M* /M on the OPEP contribution. Figure ^ 
exhibits clearly the damping due to decreasing M* /M. We stress that the mech- 
anism of damping is generic to any relativistic treatment using derivative-coupled 
pion and not particular to either Ref . |^ or the use of Ref . |^ . 
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Fig. 2. Plots of {NM I Kr I NM) with the parameters of Ref. ||]. In the G-matrix 
calculations S is replaced with aS. The plots are for a = 0., 0.5 and 1.0. The last one is 
the result of a D-B self-consistent calculation. The other two are not self consistent. 

We want to be careful that the present work not be interpreted as providing sup- 
port for MFT. Results of calculations of E/A using the same interaction, namely, 
that of Ref. [^, in both D-B and MFT treatments are shown in Fig. ||. We see 
that the results are distinctly different. Such differences are found in the results 
for scalar and vector fields in the two treatments. Results obtained upon using the 
intearction of Ref. Q are listed below. 

Sd^b = -306 (p/po)° '' MeV, 
Vd-b = 233(p/po)"'^MeV, 
Smft = -358 (p/po)°-'' MeV, 
Vmft = 295(p/po)MeV 

Undoubtedly, if one releases oneself from the constraint of fitting NN data and 
freely chooses the NN interaction one can obtain proper binding and saturation of 
nuclear matter with a MFT calculation. 

The results presented in this talk are the first explicit calculations showing that 
in the relativistic treatment the tensor force contributions generated by derivative- 
coupled OPEPHare reduced in size in nuclear matter. Because of this, in complete 




Pauli p tensor force is also damped. 
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Fig. 3. Plots of E/A from a Dirac-Brueckner and a MFT calculation with the same 
quasipotential. 



Table 2. Amorim-Tjon and Bonn C parameters for derivative-coupling and standard 
pseudoscalar coupling parameters. The value of rricr is not specified. 



Coupling 


Source 


Isospin 


ma- 


9a 












in MeV 




in MeV 




Derivative 


Amorim-Tjon |B| 




570 


7.6 


138 


14.2 




Bonn C § 


1 


550 


8.6 


138 


14.2 









720 


17.6 


138 


14.2 


Pseudoscalar 






7 


? 


138 


14.2 



contrast to the nonrelativistic situation, they cease to play an essential role in the 
saturation mechanism. The reduction of the tensor force contributions is princi- 
pally due to the relativistic M* /M effect. But even the reduced role of OPEP is 
not negligible in the actual saturation properties of nuclear matter. As noted, it 
contributes — 20MeV to E/A. The dominant mechanism of saturation of nuclear 
matter is basically very different in the two approaches. In the nonrelativistic ap- 
proach it is the density-dependent reduction due to Pauli blocking of the attraction 
from tensor force, while in the relativistic approach it is the reduction of the rate 
of growth with increasing p of the attraction from the scalar field relative to the 
growth of repulsion from the vector field. 

Finally, let us discuss the issue of derivative versus pseudoscalar ttN coupling. 
If the former is the correct coupling for nuclear effective lagrangian then nonrela- 
tivistic treatment of nuclear matter appears not to be valid. On the other hand, if 
a pseudoscalar coupling lagrangian can be found which gives satisfactory results for 
nuclear matter the nonrelativistic treatment may be valid. Unfortunately there are 
no published results for nuclear matter with a pseudoscalar coupling lagrangian. 
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Needless to say, before doing any nuclear matter calculation the parameters must 
be fixed by fitting NN data. 

It is useful to remind ourselves that we have been dealing with quasipotentials. 
These are still constrained to be chiral invariant. The interaction lagrangians with 
which one could reproduce the quasipotentials via tree graphs with form factors 
are listed below. 

^-0757^fV' • 5^7? + ga'ipcrip. (8) 
g-rripiff + i"t5T ■ T:]'ip. (9) 

In the derivative-coupling lagrangian the nucleon fields are unaffected by chi- 
ral transformation and the scalar field cr is a chiral singlet. In the pseudoscalar 
coupling lagrangian the nucleon fields, ip and 'ip, belong to chiral (1/2,0), (0, f/2) 
representations, while a and 7? form chiral (1/2,1/2) representations. In Table g 
we list the parameters for two derivative-coupling lagrangians, namely, Amorim- 
Tjon 1^ and Bonn C Q. The last line gives the standard pseudoscalar coupling 
lagrangian parameters. Notice that both the mass and the coupling constant of 
the a meson have been left unspecified. The reason is that, quite unlike OPEP, 
there will be considerable modification of the one-u exchange potential as one goes 
from the form given by the original lagrangian to the quasipotential. This happens 
through two distinct mechanisms. First, the a couples to the pion clouds of each 
of the pai! r of interacting nucleons. Second, a pair of interacting pions are ex- 
changed between two nucleons. The interaction must be isovector in the two-pion 
t-channne|^. The prospect of a pseudoscalar lagrangian succeeding in the nuclear 
matter problem is not very good. The undamped tensor force will contribute an 
additional ~ — 15MeV. To compensate this the effective a nucleon coupling in the 
quasipotential must decrease. It is difficult to see how such a reduction can come 
about. Still, the only recourse is to actually carry out the program of study with a 
pseudoscalar lagrangian before we come to a definitive conclusion about the future 
nonrelativistic treatment of nuclear matter. 

This research is supported by the U.S. Dept. of Energy under grant no. DE- 
FG02-93ER-40762. 



^derivative coupling 
'-ttN 

fli, coupling 
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